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The cerium a-j phase transition is characterized by means of a many-body Jastrow-correlated 
wave function, which minimizes the variational energy of the hrst-principles scalar-relativistic Hamil¬ 
tonian, and includes correlation effects in a non-perturbative way. Our variational ansatz accurately 
reproduces the structural properties of the two phases, and proves that even at temperature T = OK 
the system undergoes a first order transition, with ab-initio parameters which are seamlessly con¬ 
nected to the ones measured by experiment at hnite T. We show that the transition is related to 
a complex rearrangement of the electronic structure, with key role played by the p-f hybridization. 
The underlying mechanism unveiled by this work can hold in many Ce-bearing compounds, and 
more generally in other f-electron systems. 

PACS numbers: 71.15.-m, 71.20.Eh, 71.27.+a, 02.70.Ss 


Understanding the anomalous behavior of cerium, the 
prototypical /-electron system, is one of the main chal¬ 
lenges in condensed matter physics. The 4/ electrons are 
strongly localized and their on-site Coulomb repulsion is 
large compared to bandwidth. Among all lanthanides, 
cerium is particularly fascinating, due to the strong hy¬ 
bridization with the 6s6pbd bands, all present at the 
Fermi level. The origin of the cerium volume collapse 
along the isostructural a-j transition has been a puzzle 
since its discovery in 1927[T]. A microscopic comprehen¬ 
sive description of the transition is still lacking, because a 
direct comparison with the measured structural proper¬ 
ties requires an accuracy below 10 meV. This challenges 
any ab-initio method, particularly in a regime of strong 
correlation. Model calculations have been performed in 
the Mott [2], Kondo[3l[4] and dynamical mean field the¬ 
ory (DMFT) [5H3 frameworks, with input parameters 
either chosen ad-hoc or derived from first-principles den¬ 
sity functional theory (DFT) and cRPA calculations pT]. 
Fully first-principle electronic structure schemes, such as 
DFT[n] or Gwng, grasp some features of the a and 7 
phases, but the quantitative agreement with experiment 
is generally quite poor. 

Experimentally, pure cerium undergoes the a-j tran¬ 
sition always at finite temperature T. Recently, very ac¬ 
curate X-ray diffraction measurements undoubtedly con¬ 
firmed the first-order Fm3m isostructural character of the 
transition [ 13 ]. The first-order line extrapolates to zero-T 
at negative pressures. Nevertheless, the T=0K determi¬ 
nation of its phase diagram is extremely important as 
it can shed light on the underlying electronic structure 
mechanism of the transition, and clarify some critical 
points still under debate. For instance, some experiments 
with cerium alloys seem to find a critical low-T end-point 
on the (a -7 phase boundary [14], where the effect of alloy¬ 


ing is expected to provide a negative chemical pressure on 
the cerium sites. However, it has also been proven that 
the end-point of the critical line can be tuned down to 
zero T by changing the bulk modulus through alloying, 
thus opening the way of new low-T scenarios, like super¬ 
conducting and non-Fermi liquid fluctuations [15]. The 
presence of a low-T end-point is obviously material de¬ 
pendent and it is therefore possible that cerium allows 
instead a genuine /-electron driven OK quantum phase 
transition in the negative pressure side of its phase dia¬ 
gram. 

In this paper, we present a detailed analysis of the elec¬ 
tronic structure modification across the volume collapse, 
studied from first-principles, by means of an explicitly 
correlated many-body wave function and accurate OK 
quantum Monte Carlo (QMC) techniques. Remarkably, 
we have been able to stabilize two distinct coexisting so¬ 
lutions, a and 7 , with the full set of structural parameters 
across the transition seamlessly connected to the exper¬ 
imental values at finite T. We prove that the transition 
results from a subtle competition between local Coulomb 
repulsion and bandwidth, the latter determined mainly 
by the aig and Uu atomic orbitals. The key role is played 
by the p-f hybridization, set by the octahedral crystal 
field, which allows the tiu orbital to breath between the 
two phases. In the 7 phase, the chemical bond has weaker 
aig and stronger tiu channels, due to more extended tiu 
orbitals, if compared to the a phase at the same volume. 
This weakens the bond strength while it reduces the on¬ 
site Coulomb repulsion, resulting in a stabilization of the 
7 phase at larger volumes. 

In our approach, the two phases are described by 
a paramagnetic Jastrow-correlated Slater determinant 
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(JSD) wave function sampled by QMC techniques: 

^jSD(Rei) = exp[-J(Rei)] det[0i(rj)], (1) 

where 1 < i, j < N, Rei = {ri ..., r^v} is the many-body 
A/'-electron configuration, and the determinant is factor¬ 
ized in two spin components t and J., since the molecular 
orbitals (pi have a definite spin projection along z. Both 
J and p are analytic functions with parameters that min¬ 
imize the energy of the scalar-relativistic first-principles 
Hamiltonian (see [16] for details). The full Coulomb 
electron-ion interaction is replaced by a scalar-relativistic 
Hartree-Fock energy consistent pseudopotential [17] with 
5<s^5p^6<s^5d^4/^ atomic reference configuration, which 
includes semi-core states. 

The Jastrow factor takes into account strong local cor¬ 
relations as well as intersite correlations, and thoroughly 
modifies the DFT generated Slater determinant. We fully 
optimized the JSD wave function in a 32-atom cubic su¬ 
percell with periodic boundary conditions, which yields 
structural parameters close to the thermodynamic limit 
(see Supplemental Material in m)- By determining 
the variational energy as a function of the unit cell vol¬ 
ume, we evaluated the equation of states at the varia¬ 
tional Monte Carlo level (VMC), as reported in Fig.j^a). 
The fee equilibrium volume per atom of the a phase 
turns out to be 27.dA^, in a greatly better agreement 
with the experimental value (28.52A^[T8]) than LDA or 
GGA DFT calculations. To further improve the elec¬ 
tronic structure, we used the lattice regularized diffusion 
Monte Garlo (LRDMC) method [T9ll2Q] . In the LRDMG, 
the starting point is our best VMC wave function given 
byQ, that is projected to the ground state with the ap¬ 
proximation of the fixed nodes, pinned to the ones of the 
VMC wave function to cope with the sign problem arising 
in the imaginary time projection. The LRDMC equation 
of states plotted in Fig. [^b) yields an equilibrium vol¬ 
ume of 28. 4 A^, in very good agreement with experiment, 
while the bulk modulus B is overestimated (see Tab.[^. 



VMC (T=0K) 

LRDMC (T=0K) 

exp 

(A^) 

27.4 ± 0.1 

28.4 ± 0.2 

28.52[TH] 

Wq (A^) 

30.8 ± 0.2 

32.3 ± 0.3 

34.35[1I] 

B“ (GPa) 

48 ± 1 

50 ± 3 

35 [22] 

B'>' (GPa) 

38 ± 1 

45 ± 3 

21-24(231211 

K.in (A") 

28.0 ± 0.2 

28.5 ± 0.3 

28.2[I3] 

K.ax (A^) 

31.3 ± 0.3 

32.7 ± 0.4 

32.8(13] 

5V{%) 

11.7 ± 0.6 

13.8 ± 1.1 

15.1(13] 

Pt (GPa) 

-0.63 ± 0.29 

-0.45 ± 0.53 

0.7(13] 

AU (meV) 

13 ± 1 

12 ± 3 

25(13] 


TABLE I: Structural and phase transition parameters for 
a and 7 phases obtained by VMC and LRDMC, compared 
with the experiment. The 0-7 phase transition parameters 
are taken from Ref. [13] at T = 334K. A further detailed com¬ 
parison with alternative ab-initio methods is reported in the 
Supplemental Material. m 


By starting from the optimal JSD wave function for 
the a phase, we performed VMC energy minimizations 
at much larger volumes (> 40A), where we stabilized a 
second paramagnetic solution, lower in energy than the a 
phase. This second solution holds out even at smaller vol¬ 
umes, although at higher energies. Our computer simu¬ 
lations then reproduced what is seen in experiment, with 
a clear hysteresis between the two states as a function 
of volume (see Fig. [^c)). Further analysis, based on the 
Maxwell common tangent construction and on the ge¬ 
ometry parameters, confirmed that this solution is fully 
compatible with the sought 7 phase. Its LRDMC equi¬ 
librium volume is 32.3A^ (Tab. [^, representing a vol¬ 
ume collapse of about 13%; its bulk modulus is softened 
with respect to the a phase (as seen in experiments); the 
(negative) transition pressure pt is compatible with the 
extrapolated transition line to the negative side of the ex¬ 
perimental p-T diagram; and the lower Vmin and upper 
Vmax critical volumes are within the experimental range 
of phase coexistence (Fig. [^c) and Tab.|I|. 

Once the macroscopic parameters are determined, our 
theoretical approach is qualified to provide the micro¬ 
scopic physical origin of the volume collapse transition. 
By QMC methods it is actually possible to directly ac¬ 
cess spin and charge fluctuations, through the measure 
of the spin-spin and charge-charge correlation functions. 
Here, we define the charge and spin operators on a cerium 
site as hi = fv(^R-)dr + q(r)'*/)^(r)) and 

di = (q(r) 7 (r) - q(r)?/)^(r)), where the 

fermionic field' 0 j(r) (' 0 ^(r)) creates (annihilates) an elec¬ 
tron of spin a at the position r, and the integral is done 
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FIG. 1: Panels (a) and (b): Equation of states of the a 
and 7 phases obtained by VMC and LRDMC calculations, 
respectively. The black dashed straight line is the Maxwell 
construction with the corresponding calculated 0-7 transition 
pressure. Panel (c): Clapeyron diagram obtained at OK in 
quantum Monte Carlo, compared to the experimental phase 
diagram by Decremps et al. m at finite T. Remarkably, the 
upper and lower critical volumes are within the experimental 
range of the coexistence region. 
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over a sphere of radius R = 2.5 a.u. around the nu¬ 
cleus Hi. At the volume V=31.73 A^, which falls into 
the experimental coexistence region, this integration ra¬ 
dius gives (hi) ~ 9 electrons per site in both phases, 
mainly coming from the semi-core and the 4/ 

states, which are the most localized among the valence 
electrons described by our pseudopotential. The on-site 
charge fluctuations {fiifii) — (hi){hi) computed by VMC 
are 1.32(2) and 1.35(1) for the 7 and a phase, respec¬ 
tively. There is no sizable difference between the two 
phases. LRDMC does not change this picture. More¬ 
over, the Jastrow parameters which control the charge- 
charge correlations do not change significantly between 
the two phases, in accordance with the charge-charge cor¬ 
relation function results. Therefore, no suppression of 
double occupancies occurs in the 7 phase, signaling that 
the Mott scenario of the a-j transition should be defi¬ 
nitely discarded. This is an important conclusion, con¬ 
sidering that the Mott transition has been proposed as 
a valid interpretation of the volume collapse until very 
recently [ 2 j |TTJ [25] . 

In the spin sector {hi) = 0, because the a and 7 wave 
functions do not break the spin symmetry, as both states 
are paramagnetic by construction. From the experimen¬ 
tal point of view, the cerium fee lattice undergoes the vol¬ 
ume collapse between two paramagnetic states at finite 
temperature. However, the a and 7 phases feature a very 
different magnetic susceptibiiity, the former being Pauii- 
iike, the fatter of Curie-Weiss type. Eariy caicuiations 
based on the Kondo modei [3] and iater LDA+DMFT 
studies [IHH] expiained this difference in terms of Kondo 
iocai moment formations in the 7 phase, whiie the ef¬ 
fective Kondo temperature of the DMFT impurity prob- 
iem is much iarger in the a phase, ieaving it in the fuiiy 
screened singiet state. Thus, the spin response charac¬ 
terizes the two paramagnetic phases at finite tempera¬ 
ture. At OK the on-site spin-spin correiation functions 
{hihi) — {hi) {hi) computed by VMC yieid 0.5614(3) and 
0.5861(4) for the 7 and a phase, respectiveiy. Moreover, 
the spin-spin correiations have a very short-range. As 
in the charge sector, there is no significant difference be¬ 
tween the two phases. This resuit can be understood in 
terms of Kondo physics. T= OK is iower than any finite 
Kondo temperature, so that both phases are in the fuiiy 
screened regime [26]. 

The most striking difference between the a and 7 so- 
intions is in the eiectron density {p{x^y^z))^ computed 
aiways the the same voiume V=31.73 A^. The xy con¬ 
tour piot of pa— P'y is presented in Fig.j^a) for 2 : = 0, i.e. 
at the piane containing the centrai atom and 4 nearest 
neighbors. This shows positive (red) and negative (biue) 
iobes of atomic character. The Jastrow factor cannot 
expiain this difference on its own, as we have seen that 
its variation between the two phases is weak. Instead, 
the difference shouid come from the Siater determinant 
= det[(/)], once it is combined and optimized to¬ 


gether with the Jastrow factor in 0 . This is indeed con¬ 
firmed by p^ — where p^^ = with J 

dropped (set to 0) and frozen, plotted in Fig. |^b). 
The charge density difference carried by the Slater de¬ 
terminant shows a similar pattern than the full many- 
body JSD density. This is a strong signature that the 
main difference between the a and 7 phases at OK comes 
from a statie rearrangement of the electronic structure, 
driven by the dynamie electron correlation, which - in 
our approach - is coded in the Jastrow factor. It is worth 
pointing out here that both phases share almost the same 
radial charge density (up to a 0.5% difference) around 
the nuclei. The main variation is in its angular distribu¬ 
tion, suggesting that the transition must be understood 
in terms of an electron rearrangement at the atomic level, 
which will consequently affect the chemical bond in the 
solid. 

Ji ;7 
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(c) ' (d) 

FIG. 2: xy contour plot at 2 : = 0 (plane containing 

the central atom). Panel (a): JSD-VMC density difference 
Pa — P-f] Panel (b): p^ — p^^ density difference coming 
from the determinantal part only of the JSD-VMC wave 
function. Panels (c), (d): Spf^ = X)j {~ ? with 

the projected electron density, 
where i = {s + d, p + /(tin) + /(t2n)} are the atomic or¬ 
bital symmetries for the (c), and (d) panels, respectively. The 
density values are expressed in A“^. The JSD-VMC density 
values are twice smaller than the color-code scale printed in 
the key. The unit cell volume is 31.73A^. The location of Ce 
atoms are indicated by gray dots. The nearest neighbors on 
the plane of the central atom are at the square corners. 

In order to analyze this hypothesis, we consider the 
density matrix D^^^^{r^r') = 

jected over a single cerium atom. This is obtained by ex¬ 
panding the molecular orbitals ' 0 i(r) on an atomic basis 
set and considering in (r) only the components refer- 
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ring to the chosen atom. We then determine the “best” 
atomic orbitals representing the projected den¬ 
sity matrix by optimally re¬ 

duced space, namely in terms of only k « N atomic nat¬ 
ural orbitals (ANOs) centered on the reference atom and 
corresponding auxiliary molecular orbitals span¬ 

ning all the cell. This can be achieved by a standard 
Schmidt decomposition, through a minimization of the 
Euclidean distance between the truncated and the pro¬ 
jected density matrix. The resulting eigenvalues Af are 

such that ^ Yli=i^h related to 

the ANOs occupation and their chemical reactivity (see 
Supplemental Material in m)- 



FIG. 3: Left panel: First 13 ANO eigenvalues of the density 
matrix Schmidt decomposition for the a and 7 phases at E = 
31.73 A^. Their symmetry is reported in the x-axis. The 
vertical red line represents the filling ( 6 ) of a non-interacting 
closed-shell pseudoatom in the octahedral field. The variation 
between the two phases is remarkable for the aig (atomic 
6 s) 7-th eigenvalue, which is more resonating with the ones 
below in the a phase. Right panel: natural orbital localization 
measured by the spread Q = (T|r^|T) — |(T|r|T)|^. Note that 
the largest difference comes from the tiu orbitals. 

Any local atomic variation due to a change in the chem¬ 
ical bond or crystal field is detected by this approach, as 
it takes into account the embedding of the atom in its 
environment. In Fig. we plot the ANOs A? eigenvalues 
and their spread, for the same volume as in Fig. The 
first 6 ANOs would be perfectly occupied in case of non¬ 
interacting closed-shell pseudoatoms. In particular, the 
first 4 are the semi-core states. 

Two striking features are apparent. Firstly, between 
the a and 7 ANOs, there is the variation of the 7-th 
atomic orbital weight. In the a phase the 7-th ANO, of 
aig symmetry, has almost the same weight as the 5-th 
and 6 -th ANOs, of Cg symmetry arising from 
and 5 dx 2 -y 2 orbitals, degenerate in the octahedral crys¬ 
tal field. As the isolated atomic ground state is in the 
singlet Gs^bd^Af^ configuration (with the s shell full 
and inert), it is clear that the s-to-d atomic promotion 


is crucial to explain the chemical bond in the a phase, 
with the 2 eg (5d) and the aig ( 6 s) orbitals cooperating 
to set its strength. The cooperative action of s and d 
orbitals has been highlighted also in the formation of the 
Ce dimer [271128]. 

On the other hand, in the 7 phase, the aig ( 6 s) ANO 
weight is almost an order of magnitude smaller. It means 
that the s-character of the outer-shell region is weaker 
in 7 than in a, as shown also by the density variation 
projected onto the s-atomic orbitals, plotted 
in Fig.jl^c). 

In both a and 7 phases, above the 7-th aig ANO, there 
is a series of orbitals with small eigenvalues but com¬ 
peting each other in magnitude. The remaining atomic 
f-electron shares a mixed character, with non-negligible 
^ 2 u and tiu occupations. 

The second important difference between a and 7 is 
the localization of the three degenerate tiu orbitals. The 
octahedral crystal field makes the tiu orbitals strongly 
hybridized between the p and / atomic symmetries. It 
turns out that the spread of the tiu is twice larger in 
the 7 phase, as reported in Fig. [^b). This is consistent 
with Fig. I^d), where the p + f{tlu) projected density 
is spread over a wider range in the 7 phase. The larger 
tiu extension reduces the strong local Coulomb repulsion 
and increases the overlap with its neighbors and thus its 
bonding character. 

To summarize, the sizable difference between the a and 
7 chemical bond character results from a reduction of the 
ai^ weight together with an increase of the tiu overlap. 
The breathing of the tiu orbitals takes place through the 
hybridization between the p and / states, coupled by 
the octahedral crystal field. The chemical bond in the 7 
phase is weaker (and so the equilibrium volume is larger) 
as the tiu-based bond is less strong than the s-aig one. 
On the other hand, the system gains energy by reducing 
the on-site Coulomb repulsion through more extended tiu 
orbitals. 

In conclusion, the volume collapse transition can be 
understood at OK as a conventional first order transition 
of electronic origin. The two phases are well described by 
the zero-T equation of state, while their relative stability 
is provided by tiny entropic effects [13]. The underly¬ 
ing mechanism of the volume collapse should survive by 
the addition of the spin-orbit coupling [29], not present 
in our calculations, as in cerium it is much weaker than 
the local Coulomb repulsion, although competing to the 
crystal field splitting [30]. Our picture disproves the va¬ 
lidity of the Mott model, and puts cerium in a quantum 
phase transition regime. Our detailed predictions on the 
interplay between valence and localized orbitals can be 
experimentally tested by X-ray electron spectroscopy at 
L edges, to probe the s and p states. The electronic 
phase transition mechanism detailed in this work can be 
applied to cerium alloys, and more generally to a wider 
class of /-electron systems. 
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